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Abstract 

The steady state of ions diffusion in polymer electrolytes at arbitrary applied voltage is analyzed 
in the framework of the Nernst-Planck-Poisson equation (NPP). The exact solution of the set of 
equations is found without the assumption of low ions concentration. The solution is independent of 
the kinetic properties of the system. At constant voltage there is a master curve for concentration in 
terms of the initial concentration. Enhancing the voltage causes an increase of the ion concentration 
gradient and consequently the applicability of the NPP is violated for high voltages. The analytical 
finding is estimated by using experimental data from recent measurements (P. Kohn et al Phys. 
Rev. Left. 99, 086104 (2007)). As the result we find an upper bound for the vahdity of the NPP. 
Above this voltage higher order gradient terms become relevant. 
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I. INTRODUCTION 



Polymer electrolytes characterized as amorphous structure are widely used for solid state 
ionics devices such as lithium ion batteries fuel cells, electrochemical displays etc. The 
material is a composite one consisting of dissociated ions which are dissolved in a polymer 
or glassy matrix. Very recently the motions of proteins in densely grafted layers of polyelec- 
trolyte brushed has been studied 21] or the dynamics of polyelectrolyte transport through a 
protein channel as a function of applied voltage [3(. The understanding of the ionic transport 
mechanism in polymer materials is of basic interest and has attracted attentions over several 
decades, for a recent review on the history, the applications of ionic transport and structure 
as well as dynamics see jj]. A number of different theoretical methods has been proposed 
for polyelectrolyte systems . The activities concerning the theoretical approach for 

the ionic transport behavior had been reviewed in 8|. 

Our analysis is based on the Nernst-Planck equations supplemented by the static Poisson 
equation. These set of equations, abbreviated as NPP, are derived for a low charge carrier 
density. The diffuse charge dynamics is characterized by a linear gradient term for the 
ionic concentration. In addition to the diffusive charge dynamics the transport process is 
driven by an external field E. For low voltages the linearized NPP can be solved by Laplace 

n n 

transforms [81 whereas for high voltages numerical solutions had been obtained in [|8|]. As 
estimated in 8| the validity of the linear approximation should be guaranteed up to 70 mV. 
A detailed analysis of the non-linear regime is still lacking. Otherwise, the mobility of ions 
has been determined recently by transient current measurements at high voltages up to the 
order of 300 V 9|. The problem arisen is whether the NPP with a linear concentration 
gradient are also valid for such high fields. Intuitively one expects that high voltages should 
give rise to higher order concentration gradients of the mobile ions. Otherwise, very low 
voltages cause only low gradients even in systems with a high concentration of ions. Thus, 
the limitation of the NPP is related to both high concentration of charge carriers as well 
high applied voltages. However our analysis offers that a high initial concentration should 
be irrelevant. Whereas all previous papers are mainly devoted to a numerical solution of 
the NPP or some modifications of those equations such as by including steric effects 10|, 
the present paper is directly addressed to the NPP, in order to check the ability of the 
equations at high voltages. Despite the model is used widely it seems to be difficult to 
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estimate the condition under which the NPP are vahd. So far, the NPP has not been solved 
exactly but instead of that the solution had been discussed in a linear regime where the 
long time behavior was found. Here we present the exact steady state solution of the full 
NPP. Especially we demonstrate that the solution within the linear regime may be spread 
up to 70 mV approximately. Otherwise the experimental measurement of the mobility is 
fulfilled in between 70 and 300 V [9]. Obviously the NPP is not appropriate to clarify the 
observations made in j^. Apparently, a high voltages requires higher order diffusive terms. 



II. THE MODEL AND THE SOLUTION 



Let us consider a solution of charged particles embedded in polymer electrolyte at tem- 
perature above the glass-transition temperature of the polymer. We study completely dis- 
sociated electrolyte placed in between infinite flat electrodes at coordinates x = —L and 
X = L. Therefore the system offers a simple one- dimensional geometry. Concentrations of 
the charged ions are described by continuum fields C±{x, t). The charged ions are subjected 
to an electric field E = — V0. The scalar potential obeys the Poisson equation 

-e^ = .e(C+-C_), (1) 

where e is the dielectric permeability, e is the elementary charge and ze is the charge of 
the ion. Because the charge is conserved the time derivative of C± is given by the spatial 
derivative of the corresponding ionic fluxes denoted by F±. Generally this flux is composed 
of two terms, namely the diffusive part originated by a small concentration gradient and the 
electric field term. This field is assumed to couple linearly to the concentration field. The 
set of evolution equations reads 

^ = -— = -— (-D^ - zeC+^ 
dt dx dx \ dx dx 

^ - -—(F ) - -— (-D^ + fizeC (2) 
dt dx dx \ dx dx J 

Here D is the diffusivity and fi denotes the mobility, for simplicity we have assumed = 

D_ = D. Further it is supposed that the Einstein relation D = fikT is fulfilled. Likewise it is 

supposed that the Faraday currents can be neglected. The boundary conditions imposed are 

the disappearance of fluxes F± = at x = ±L. This set of equations combined with Eq. ([1]) 

are called Nernst-Planck-Poisson equations (NPP). Essentially is that Pick's law of diffusion 



is fulfilled, i.e. the ionic flux is proportional to the first derivative of concentration. At 
high voltage, the gradient of concentration becomes more pronounced and the consecutive 
equations should be changed due to higher order gradient terms. To estimate the validity 
of the NPP for high voltages let us consider the steady state. This state is characterized by 
constant ionic fluxes where the constant is zero due to the boundary conditions. For that 
case we found an exact solution demonstrated further. The difference of the charge carriers 
due Eqs. ([2]) can be expressed by the potential of the electric field according to Eq. ([1]). It 
results 

— ^^^-/i^e^^C(x) = with C{x)=C+{x) + C.{x). (3) 
Using Eqs. ([2]) and ([T]) the sum of the concentration C{x) satisfies 

where r is an integration constant. Both Eqs. ([3]) and (jlj) can be combined to a single 
equation for the electric potential 
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0. (5) 



ze dx^ D dx 

Then the solution of the last equation can be applied to find the charge concentration via 

d\nC±{x) fized(f){x) , . 

dx ~^ D dx 

Actually the solution of Eq. ([5]) is given in terms of Jacobian elliptic functions [l^. To that 
aim let us substitute y{x) = in Eq. After some formal steps we end up with an 

expression for the electric field E{x) 

E{x) = -^^ = q-(3-sn{Px-Xo,q). 7 

dx fize 



12| we get the stationary con- 



Using Eq. and the properties of the Jacobian functions 
centration profiles of the charge carriers 

C±{x) = — -q- 0^ ■ [dn{(3x -Xo,q)±q- cn{(3x - Xq, q)f . (8) 
jjLz^e^ 

The quantities sn, cn and dn are the Jacobian elliptical functions, which are characterized 
by two parameters denoted as j3 and q. Due to the symmetry property E{x) = E{—x) the 
integration constant Xq should be equal to quarter of the period of Jacobian sinus Xq = K{q), 



where K{q) is the complete elhptic integral of the first kind [12!]. The parameters (3 and q 
are determined by the following conditions: 

-L 



C±{x)-dx = 2Lr] 

-L 

E(^L) = -^. (9) 

The first condition determines the total concentration t] which is given initially. Thus, rj is 
the initial particle concentration before the voltage were applied. The second condition fixes 
the applied voltage V. After performing the integrating of the first equation in Eq. and 
reorganization of the the second one we have 

E(0L,g)-<,'m{PL,g)cd(,3L,<,) ~ ^ (1 - .f) ^ -^J^ 

Here E{f]L, q) stands for the incomplete elliptic integral of the second kind and we have 
introduced 

,fn ^ - cn{f3x,q) 

Eqs. (I7j), ([HD together with Eqs. ffTU]) provide the complete solution of the stationary problem 
of the NPP. However the conditions in Eqs. ffTOj) for the parameters jS and q are too compli- 
cated to get analytical results. To proceed one can use the fact that in the first equation the 
first term is slowly changing almost linear function, whereas the second term can be roughly 
approximated by just q. This approximation leads to an parameter area which enables us 
to get reasonable values for the parameters. 

III. RESULTS 

All calculations were performed for a real system described in 9j. The experiments had been 
carried in thin films with L = 5 ■ 10~^ m. The initial concentration rj is varied from 1 • 10^^ 
to 1.3 ■ 10^^ m~^, whereas the applied voltage V was changed from 1 • 10~^ to 2.6 • 10^^ V. In 
case the parameter q of the Jacobian elliptic function is restricted to the internal < g < 1, 
the functions are real. Inserting the experimental data obtained in [o] in our solution we find 
g >> 1 and (3 << 1. In particular it results q ~ 10^ and f3 ~ 10^^. For such values of the 
parameters q and (3 the Jacobian functions become complex functions. However it appeared 



that << and QC± << ^C±. We estimate that the imiginary part is about 7 or 8 
orders of magnitude smaller than the real part. Despite of the complexness of the Jacobian 
functions, the solution found in the paer has a physical sense. 



A. Results for constant voltage 

Calculations at constant voltage V, defined in Eq. (j9]), offer that the parameter q changes 
proportional to the initial concentration rj, whereas the other parameter/5 behaves propor- 
tional to r]~^, i. e. we can set /3 = ^ and q = qr], where the tilde means a value at some 
definite concentration. Notice that for different applied voltage V the parameters q and (3 
are also different. As one can see from Eq. ([9]), the integral J^j^ ^^^^^ ■ dx is independent 
on the fixed initial concentration i] indicating the existence of a master curve for Our 
calculation confirms the presumption. In FigH] the spatial charge distribution is shown ac- 
cording to Eq. (IH]). The different curves depicted in Fig. [T] correspond to different values 
for rj. In FigE] one observes that all these data for the concentration in terms of rj lay on 
a single master curve. Further discussions concerning scaling properties of the curves one 



can in ll| . The profile of the electric field is independent of the concentration which is 
demonstrated in Fig. [3l with other words the efficiency of screening is independent of t]. 
This observation would suggest the applicability of NPP is guaranteed also for a high charge 
carrier concentration. According to Eq. ([6]) a constant electric field suggest an increasing 
of the concentration gradient proportional to the concentration itself. Therefore high gra- 
dients could lead to an inconsistency of NPP. Another feature mentioned is the so called 
charge inversion, i.e. interfacial charges attract counterions in excess of their own nominal 
cha.,e. This effec. we. found bo.h in theo.e.ical Q and expenn,en.a. QQ wo.U. fo. 



an overview and discussion see 
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But, as one can see from Figj2l in the framework of 



steady state solution of NPP this effect is absent. 



B. Results for constant concentration 



For different voltages but constant concentration there exists no master curve as it is shown 
in Fig. m For increasing applied voltage V we find oscillations of the electric field E{x). 
From here we conclude that the NPP equations are not longer adequate in describing the 
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physical situation in mind. Such oscillations appear when the period of the Jacobian sinus 
is sufficiently small. Analytically the criterion for the absence of oscillations is the condition 
f3L < K{q). In case of (3L = K{q) the half period of the Jacobian sinus fits exactly the 
capacitor between — L < x < L. To be more specific we have checked that for an arbirtrary 
but fixed initial concentration rj the inequality f3L < K{q) is only satisfied in case the applied 
voltages fulffils V < 0.22 V. Because for V > 0.23 V the inequality is violated we conclude 
that there exist a kind of critical applied voltage in between 0.22 < V < 0.23 V. If any the 
dependence of that result on the initial concentration is weak. But such a weak dependence 
should exist due to Eq. ([6]), which signalized that the concentration gradient is proportional 
to both the applied voltage V and the initial concentration t]. 

The steady state solution of NPP equations is universal and nearly independent on specific 
system properties. Thus the kinetic properties should depend on the diffusion coefficient D. 
But in the expression for the steady state according to Eq. ([8]) there appears only the ratio 
D/fi. If Einstein relation is fulfilled the solution is independent on the kinetic properties of 
the system but depends only on the temperature. In other words, different systems with the 
same ze tend to the same steady state solution, which is however reached on different ways. 
If the voltage V is sufficiently high to alter the transport properties of polymers (mechanisms 
for a change are discussed in [ITi]) but not high enough to alter its dielectric permeability, 
the steady state will be unchanged. 



IV. CONCLUSION 



By getting the steady state we found a limit of the applicability of the Poisson-Nernst- 
Planck equations for both high voltages and high concentrations. Using the experimental 
data from we could estimate that above V = 0.22 V the NPP yield unphysical results 
characterized by the appearance of unphysical oscillation of the solutions. Notice that 
the validity of the linear approximation is restricted to an applied voltage V < 0.07 
V js]. The limit of the applicability of NPP due to a high concentration of ions is 
not very sharp. Unlike a behavior at high voltage, in this case NPP do not contain 
explicit restrictions itself. Existents of the master curve for normalized concentration is 
an indication for that fact. The inclusion of higher order gradient terms leads a more 
complicated model. Otherwise the assumption that the flux is proportional to a linear gradi- 
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ent term is a very promising assumption which should be vahd for the majority of situations. 

One of us (A. G.) acknowledges support by the International Max Planck Research School 
for Science and Technology of Nanostructures in Halle, Germany and Thomas Michael for 
helping in preparing the figures. Further we thank Peter Kohn for useful discussions. 
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FIG. 1: Concentration profile C+{x) in the capacitor. Different curves correspond to different 
initial concentrations rj. The applied voltage V is assumed to be 5 • 10~^V. 
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FIG. 2: Concentration profile ^^^^^ in terms of the initial concentration rj in the capacitor. The 
applied applied voltage is y = 5 ■ lO^^V. 
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FIG. 3: Electric field profile E{x) within the capacitor. The profile is independent on the initial 
concentration rj. The voltage is assumed to be 5 • 10~^V. 
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FIG. 4: Concentration profile C+(,t) versus coordinate in the capacitor. Different curves correspond 
to different applied voltages V . The initial concentration r/ is 5 • 10^^. 
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